FERMILAB-PUB-12-574-E 



Measurement of the differential 7 + c-jet cross section and the ratio of differential 
7 + c and 7 + b cross sections in pp collisions at \/s = 1.96 TeV 

V.M. Abazov, 32 B. Abbott, 69 B.S. Acharya, 26 M. Adams, 47 T. Adams, 45 G.D. Alcxccv, 32 G. Alkhazov, 36 

A. Altom 1 , 58 A. Askew, 45 S. Atkins, 56 K. Augsten, 7 C. Avila, 5 F. Badaud, 10 L. Bagby, 46 B. Baldin, 46 
D.V. Bandurin, 45 S. Banerjee, 26 E. Barberis, 57 P. Baringer, 54 J.F. Bartlett, 46 N. Bartosik, 39 U. Bassler, 15 
V. Baztcrra, 47 A. Bean, 54 M. Begalli, 2 L. Bellantoni, 46 S.B. Bcri, 24 G. Bcrnardi, 14 R. Bcrnhard, 19 I. Bertram, 40 
M. Bcsancon, 15 R. Bcusclinck, 41 P.C Bhat, 46 S. Bhatia, 60 V. Bhatnagar, 24 G. Blazey, 48 S. Blessing, 45 K. Bloom, 61 
A. Bochnlcin, 46 D. Boline, 66 E.E. Boos, 34 G. Borissov, 40 A. Brandt, 72 O. Brandt, 20 R. Brock, 59 A. Bross, 46 

D. Brown, 14 J. Brown, 14 X.B. Bu, 46 M. Buehlcr, 46 V. Bucschcr, 21 V. Bunichcv, 34 S. Burdin 6 , 40 CP. Buszello, 38 

E. Camacho-Pcrcz, 29 B.C.K. Casey, 46 H. Castilla-Valdcz, 29 S. Caughron, 59 S. Chakrabarti, 66 D. Chakraborty, 48 

K.M. Chan, 52 A. Chandra, 74 E. Chapon, 15 G. Chen, 54 S. Chcvalicr-Thcry, 15 S.W. Cho, 28 S. Choi, 28 
B. Choudhary, 25 S. Cihangir, 46 D. Claes, 61 J. Clutter, 54 M. Cooke, 46 W.E. Cooper, 46 M. Corcoran, 74 F. Couderc, 15 
M.-C. Cousinou, 12 A. Croc, 15 D. Cutts, 71 A. Das, 43 G. Davies, 41 S.J. de Jong, 30 ' 31 E. Dc La Cruz-Burclo, 29 

F. Dcliot, 15 R. Dcmina, 65 D. Dcnisov, 46 S.P. Denisov, 35 S. Dcsai, 46 C. Deterre, 15 K. DeVaughan, 61 H.T. Dichl, 46 
M. Dicsburg, 46 P.F. Ding, 42 A. Dominguez, 61 A. Dubey, 25 L.V. Dudko, 34 D. Duggan, 62 A. Dupcrrin, 12 S. Dutt, 24 

A. Dyshkant, 48 M. Eads, 61 D. Edmunds, 59 J. Ellison, 44 V.D. Elvira, 46 Y. Enari, 14 H. Evans, 50 A. Evdokimov, 67 
V.N. Evdokimov, 35 G. Facini, 57 L. Feng, 48 T. Fcrbcl, 65 F. Fiedler, 21 F. Filthaut, 30 - 31 W. Fisher, 59 H.E. Fisk, 46 
M. Fortncr, 48 H. Fox, 40 S. Fucss, 46 A. Garcia-Bellido, 65 J.A. Garcfa-Gonzalcz, 29 G.A. Garcfa-Gucrra c , 29 
V. Gavrilov, 33 P. Gay, 10 W. Ceng, 12 ' 59 D. Gcrbaudo, 63 C.E. Gcrbcr, 47 Y. Gcrshtcin, 62 G. Ginthcr, 46 ' 65 
G. Golovanov, 32 A. Goussiou, 76 P.D. Grannis, 66 S. Grcdcr, 16 H. Greenlee, 46 G. Grcnicr, 17 Ph. Gris, 10 J.-F. Grivaz, 13 
A. Grohsjean d , 15 S. Griinendahl, 46 M.W. Griincwald, 27 T. Guillemin, 13 G. Gutierrez, 46 P. Gutierrez, 69 J. Haley, 57 
L. Han, 4 K. Harder, 42 A. Harcl, 65 J.M. Hauptman, 53 J. Hays, 41 T. Head, 42 T. Hcbbeker, 18 D. Hedin, 48 

H. Hegab, 70 A.P. Hcinson, 44 U. Heintz, 71 C. Hcnscl, 20 I. Hercdia-De La Cruz, 29 K. Hcrncr, 58 G. Hesketh^, 42 
M.D. Hildrcth, 52 R. Hirosky, 75 T. Hoang, 45 J.D. Hobbs, 66 B. Hocnciscn, 9 J. Hogan, 74 M. Hohlfeld, 21 

I. Howley, 72 Z. Hubacek, 7 - 15 V. Hynek, 7 I. Iashvili, 64 Y. Ilchenko, 73 R. Illingworth, 46 A.S. Ito, 46 S. Jabeen, 71 
M. Jaffre, 13 A. Jayasinghc, 69 M.S. Jeong, 28 R. Jcsik, 41 P. Jiang, 4 K. Johns, 43 E. Johnson, 59 M. Johnson, 46 

A. Jonckheere, 46 P. Jonsson, 41 J. Joshi, 44 A.W. Jung, 46 A. Juste, 37 E. Kajfasz, 12 D. Karmanov, 34 P.A. Kasper, 46 
I. Katsanos, 61 R. Kchoe, 73 S. Kcrmiche, 12 N. Khalatyan, 46 A. Khanov, 70 A. Kharchilava, 64 Y.N. Kharzhecv, 32 
I. Kiselevich, 33 J.M. Kohli, 24 A.V. Kozclov, 35 J. Kraus, 60 A. Kumar, 64 A. Kupco, 8 T. Kurca, 17 V.A. Kuzmin, 34 
S. Lammcrs, 50 G. Landsberg, 71 P. Lcbrun, 17 H.S. Lee, 28 S.W. Lee, 53 W.M. Lee, 46 X. Lei, 43 J. Lcllouch, 14 
D. Li, 14 H. Li, 11 L. Li, 44 Q.Z. Li, 46 J.K. Lim, 28 D. Lincoln, 46 J. Linnemann, 59 V.V. Lipacv, 35 R. Lipton, 46 
H. Liu, 73 Y. Liu, 4 A. Lobodenko, 36 M. Lokajicek, 8 R. Lopes dc Sa, 66 H.J. Lubatti, 76 R. Luna-GarciaV 9 
A.L. Lyon, 46 A.K.A. Maciel, 1 R. Madar, 19 R. Magaha-Villalba, 29 S. Malik, 61 V.L. Malyshev, 32 Y. Maravin, 55 
J. Martinez-Ortega, 29 R. McCarthy, 66 C.L. McGivern, 42 M.M. Meijer, 30 ' 31 A. Melnitchouk, 46 D. Menezes, 48 
P.G. Mercadante, 3 M. Merkin, 34 A. Meyer, 18 J. Meyer, 20 F. Miconi, 16 N.K. Mondal, 26 M. Mulhearn, 75 E. Nagy, 12 
M. Naimuddin, 25 M. Narain, 71 R. Nayyar, 43 H.A. Neal, 58 J.P. Negret, 5 P. Neustroev, 36 H.T. Nguyen, 75 
T. Nunncmann, 22 J. Orduna, 74 N. Osman, 12 J. Osta, 52 M. Padilla, 44 A. Pal, 72 N. Parashar, 51 V. Parihar, 71 
S.K. Park, 28 R. Partridge 6 , 71 N. Parua, 50 A. Patwa, 67 B. Penning, 46 M. Perfilov, 34 Y. Peters, 20 K. Pctridis, 42 

G. Petrillo, 65 P. Petroff, 13 M.-A. Pleier, 67 P.L.M. Podesta-Lerma' 1 , 29 V.M. Podstavkov, 46 A.V. Popov, 35 
M. Prewitt, 74 D. Price, 50 N. Prokopenko, 35 J. Qian, 58 A. Quadt, 20 B. Quinn, 60 M.S. Rangel, 1 K. Ranjan, 25 
P.N. Ratoff, 40 I. Razumov, 35 P. Renkel, 73 I. Ripp-Baudot, 16 F. Rizatdinova, 70 M. Rominsky, 46 A. Ross, 40 
C. Royon, 15 P. Rubinov, 46 R. Ruchti, 52 G. Sajot, 11 P. Salcido, 48 A. Sanchez-Hernandez, 29 M.P. Sanders, 22 
A.S. Santos*, 1 G. Savage, 46 L. Sawyer, 56 T. Scanlon, 41 R.D. Schamberger, 66 Y. Scheglov, 36 H. Schellman, 49 
C. Schwanenberger, 42 R. Schwienhorst, 59 J. Sekaric, 54 H. Severini, 69 E. Shabalina, 20 V. Shary, 15 S. Shaw, 59 
A. A. Shchukin, 35 R.K. Shivpuri, 25 V. Simak, 7 P. Skubic, 69 P. Slattery, 65 D. Smirnov, 52 K.J. Smith, 64 G.R. Snow, 61 
J. Snow, 68 S. Snyder, 67 S. Soldncr-Rcmbold, 42 L. Sonncnschcin, 18 K. Soustruznik, 6 J. Stark, 11 D.A. Stoyanova, 35 
M. Strauss, 69 L. Suter, 42 P. Svoisky, 69 M. Titov, 15 V.V. Tokmcnin, 32 Y.-T. Tsai, 65 K. Tschann-Grimm, 66 
D. Tsybychev, 66 B. Tuchming, 15 C. Tully, 63 L. Uvarov, 36 S. Uvarov, 36 S. Uzunyan, 48 R. Van Kooten, 50 



2 



W.M. van Lccuwen, 30 N. Varelas, 47 E.W. Varnes, 43 LA. Vasilyev, 35 P. Vcrdier, 17 A.Y. Verkheev, 32 
L.S. Vcrtogradov, 32 M. Vcrzocchi, 46 M. Vesterinen, 42 D. Vilanova, 15 P. Vokac, 7 H.D. Wahl, 45 M.H.L.S. Wang 46 
J. Warchol, 52 G. Watts, 76 M. Wayne, 52 J. Weichcrt, 21 L. Welty-Rieger, 49 A. White, 72 D. Wicke, 23 
M.R.J. Williams, 40 G.W. Wilson, 54 M. Wobisch, 56 D.R. Wood, 57 T.R. Wyatt, 42 Y. Xie, 46 R. Yamada, 46 
S. Yang, 4 T. Yasuda, 46 Y.A. Yatsuncnko, 32 W. Ye, 66 Z. Yc, 46 H. Yin, 46 K. Yip, 67 S.W. Youn, 46 J.M. Yu, 58 
J. Zennamo, 64 T. Zhao, 76 T.G. Zhao, 42 B. Zhou, 58 J. Zhu, 58 M. Zielinski, 65 D. Zieminska, 50 and L. Zivkovic 14 

(The DO Collaboration*) 

1 LAFEX, Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, Brazil 
2 Universidade do Estado do Rio de Janeiro, Rio de Janeiro, Brazil 
3 Universidade Federal do ABC, Santo Andre, Brazil 
4 University of Science and Technology of China, Hefei, People's Republic of China 
J Universidad de los Andes, Bogota, Colombia 
''Charles University, Faculty of Mathematics and Physics, 
Center for Particle Physics, Prague, Czech Republic 
7 Czech Technical University in Prague, Prague, Czech Republic 
8 Center for Particle Physics, Institute of Physics, 
Academy of Sciences of the Czech Republic, Prague, Czech Republic 

9 Universidad San Francisco de Quito, Quito, Ecuador 
10 LPC, Universite Blaise Pascal, CNRS/IN2P3, Clermont, France 
n LPSC, Universite Joseph Fourier Grenoble 1, CNRS/IN2P3, 
Institut National Polytechnique de Grenoble, Grenoble, France 
12 CPPM, Aix-Marseille Universite, CNRS/IN2P3, Marseille, France 
13 LAL, Universite Paris-Sud, CNRS/IN2P3, Orsay, France 
H LPNHE, Univer sites Paris VI and VII, CNRS/IN2P3, Paris, France 
15 CEA, Irfu, SPP, Saclay, France 
16 IPHC, Universite de Strasbourg, CNRS/IN2P3, Strasbourg, France 
17 IPNL, Universite Lyon 1, CNRS/IN2P3, Villeurbanne, France and Universite de Lyon, Lyon, France 
18 III. Physikalisches Institut A, RWTH Aachen University, Aachen, Germany 
19 Physikalisches Institut, Universitdt Freiburg, Freiburg, Germany 
2(1 II. Physikalisches Institut, Georg- August- Universitdt Gottingen, Gottingen, Germany 
21 Institut fiir Physik, Universitdt Mainz, Mainz, Germany 
22 Ludwig-Maximilians-Universitat Miinchen, Miinchen, Germany 
23 Fachbereich Physik, Bergische Universitdt Wuppertal, Wuppertal, Germany 
2Jt Panjab University, Chandigarh, India 
25 Delhi University, Delhi, India 
26 Tata Institute of Fundamental Research, Mumbai, India 
2 University College Dublin, Dublin, Ireland 
28 Korea Detector Laboratory, Korea University, Seoul, Korea 
' 29 CINVESTAV, Mexico City, Mexico 
30 Nikhef, Science Park, Amsterdam, the Netherlands 
31 Radboud University Nijmegen, Nijmegen, the Netherlands 
32 Joint Institute for Nuclear Research, Dubna, Russia 
33 Institute for Theoretical and Experimental Physics, Moscow, Russia 
3J> Moscow State University, Moscow, Russia 
35 Institute for High Energy Physics, Protvino, Russia 
36 Petersburg Nuclear Physics Institute, St. Petersburg, Russia 
37 Institucio Catalana de Recerca i Estudis Avangats (ICREA) and Institut de Fisica d'Altes Energies (IFAE), Barcelona, Spam 

38 Uppsala University, Uppsala, Sweden 
39 Taras Shevchenko National University of Kyiv, Kiev, Ukraine 
40 Lancaster University, Lancaster LAI J^YB, United Kingdom 
''-'imperial College London, London SW7 2AZ, United Kingdom 
42 The University of Manchester, Manchester M13 9PL, United Kingdom 

^ s University of Arizona, Tucson, Arizona 85721, USA 
44 University of California Riverside, Riverside, California 92521, USA 

45 Florida State University, Tallahassee, Florida 32306, USA 
46 Fermi National Accelerator Laboratory, Batavia, Illinois 60510, USA 
J ' 7 University of Illinois at Chicago, Chicago, Illinois 60607, USA 
4 Northern Illinois University, DeKalb, Illinois 60115, USA 
J * 9 Northwestern University, Evanston, Illinois 60208, USA 
50 Indiana University, Bloomington, Indiana 47405, USA 

51 Purdue University Calumet, Hammond, Indiana 46323, USA 

52 University of Notre Dame, Notre Dame, Indiana 46556, USA 



3 



53 Iowa State University, Ames, Iowa 50011, USA 
5J> University of Kansas, Lawrence, Kansas 66045, USA 
55 Kansas State University, Manhattan, Kansas 66506, USA 
56 Louisiana Tech University, Ruston, Louisiana 71272, USA 
57 Northeastern University, Boston, Massachusetts 02115, USA 
58 University of Michigan, Ann Arbor, Michigan 48109, USA 
59 Michigan State University, East Lansing, Michigan 48824, USA 
60 University of Mississippi, University, Mississippi 38677, USA 
61 University of Nebraska, Lincoln, Nebraska 68588, USA 
fj2 Rutgers University, Piscataway, New Jersey 08855, USA 
63 Princeton University, Princeton, New Jersey 08544, USA 
6i State University of New York, Buffalo, New York 14260, USA 
65 University of Rochester, Rochester, New York 14627, USA 
66 State University of New York, Stony Brook, New York 11794, USA 
67 Brookhaven National Laboratory, Upton, New York 11973, USA 
68 Langston University, Langston, Oklahoma 73050, USA 
69 University of Oklahoma, Norman, Oklahoma 73019, USA 
70 Oklahoma State University, Stillwater, Oklahoma 74078, USA 
71 Brown University, Providence, Rhode Island 02912, USA 

72 University of Texas, Arlington, Texas 76019, USA 
73 Southern Methodist University, Dallas, Texas 75275, USA 

7J, Rice University, Houston, Texas 77005, USA 
5 University of Virginia, Charlottesville, Virginia 22904, USA 
University of Washington, Seattle, Washington 98195, USA 
(Dated: December 24, 2012) 

We present measurements of the differential cross section da/dpp for the associated production 
of a c-quark jet and an isolated photon with rapidity |y 7 | < 1.0 and transverse momentum 30 < 
p^, < 300 GeV. The c-quark jets are required to have \y* et \ < 1.5 and p!£. > 15 GeV. The ratio of 
differential cross sections for 7 + c to 7 + b production as a function of plj, is also presented. The 
results are based on data corresponding to an integrated luminosity of 8.7 fb _1 recorded with the DO 
detector at the Fermilab Tevatron pp Collider at yfs =1.96 TeV. The obtained results are compared 
to next-to-leading order perturbative QCD calculations using various parton distribution functions, 
to predictions based on the fcT-factorization approach, and to predictions from the SHERPA and 
pythia Monte Carlo event generators. 

PACS numbers: 13.85.Qk, 12.38.Bx, 12.38.Qk 



1 In hadron-hadron collisions high-energy photons are 23 

2 mainly produced directly in a hard parton scattering pro- 24 

3 cess. For this reason, and due to their pointlike electro- 25 

4 magnetic coupling to the quarks, they provide a clean 26 

5 probe of parton-level dynamics. Photons in association 27 

6 with a charm (c) quark are produced primarily through 28 

7 the Compton-like scattering process gc — > 7c, which 29 
s dominates up to photon transverse momenta with respect 30 

9 to the beam axis oipp w 70—80 GeV, and through quark- 31 

10 antiquark annihilation, qq — > 75 — > jee, which dominates 32 

11 at higher pp [l|. Inclusive 7 + c production may also orig- 33 

12 inate from processes like gg —¥ cc or eg eg, where the 34 

13 fragmentation of a final state c-quark or gluon produces 35 

14 a photon [l[ . Photon isolation requirements substantially 36 

15 reduce the contributions from these processes. Measure- 37 

16 ments of the 7 + c-quark jet differential cross section as a 38 

17 function of pp improve our understanding of the under- 39 
is lying production mechanism and provide useful input for 40 

19 the c-quark parton distribution functions (PDFs) of the 

20 colliding hadrons. 

21 In this Letter, we present measurements of the inclu- 

22 sive 7 + c-jet production cross sections using data col- 



lected from June 2006 to September 2011 with the DO 
detector in pp collisions at s/s = 1.96 TeV which corre- 
spond to an integrated luminosity of 8.7 ± 0.5 fb _1 0. 
The cross section is measured differentially as a func- 
tion of pip for photons within rapidities \y 1 \ < 1.0 and 
30 < p^, < 300 GeV, while the c-jet is required to 
have \y^ ot \ < 1.5 and pi^ > 15 GeV. In comparison 
to our previous 7 + c-jet measurement we now re- 
tain all events having at least one jet originating from a 
charm quark, as opposed to considering only the events 
in which the charm jet candidate is the jet with highest 
Pt- To increase the signal yield and study a trend in the 
data/theory ratio observed in Ref. [H, we have extended 
the rapidity [H region from |yj ct | <0.8 to |yj ct | <1.5 and 
combine regions with positive and negative products of 
rapidities, y 7 y Jct . In addition, an increased integrated 
luminosity by about a factor of nine allows the pip range 
to be extended to higher values. 

The data set and event selections used in our measure- 
ment are similar to those used in the recently published 
measurement of the 7 + 6-jet differential cross section [|| . 
However, because of the difficulty in discriminating c jets 
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from light jets, this measurement adopts a different strat-101 
cgy for the estimation of the c-jet fraction. Here we apply 102 
a significantly more stringent requirement for selecting™ 
heavy flavor jets (originating from c and b quarks) in or-™ 
der to suppress the rates of light jets (originating fromi 05 
light quarks or gluons) by an additional factor of 2.5 — 3.™ 
This small residual contribution of light jets is then sub-107 
tracted from the selected data events prior to performing™ 
the fit with the discriminant templates of 6-jets and c-jets™ 
to extract the c-jet fraction. Using this event selectionno 
criteria, we reproduce the results for the 7 + 6-jet crossm 
section, measure the 7 + c-jet cross section and calculatem 
the ratio a(j + c)/a("f + b) in bins of p T . Common cxpcri-™ 
mental uncertainties and dependence on the highcr-ordern 4 
corrections in theory are reduced in the ratio, allowing ans 
precise study of the relative (7(7 + c)/a(j + b) rates. u 6 

The DO detector is a general purpose detector dc-m 
scribed in detail elsewhere @. The subdetectors most™ 
relevant to this analysis are the central tracking sys-119 
tern, composed of a silicon microstrip tracker (SMT)i2o 
and a central fiber tracker (CFT) embedded in a 1.9 T121 
solenoidal magnetic field, the central preshower detector^ 
(CPS), and the calorimeter. The CPS is located immedi-123 
ately before the inner layer of the central calorimeter andi24 
is formed of approximately one radiation length of leadi25 
absorber followed by three layers of scintillating strips. 126 
The calorimeter consists of a central section (CC) withm 
coverage in pseudorapidity of |?7dot I < IT Qi and two end™ 
calorimeters (EC) extending coverage to |?7det| ~ 4.2, all™ 
housed in separate cryostats, with scintillators between™ 
the CC and EC cryostats providing sampling of develop-m 
ing showers for 1.1 < |%et| < 1-4. The electromagnetic™ 
(EM) section of the calorimeter is segmented longitudi-™ 
nally into four layers (EMi, i = 1 — 4), with transverse™ 
segmentation into cells of size A^dot x A</> = 0.1 x 0.1 0,i35 
except EM3 (near the EM shower maximum), where it™ 
is 0.05 x 0.05. The calorimeter allows for a precise mea-137 
surcmcnt of the energy and direction of electrons and™ 
photons, providing an energy resolution of approximately 139 
4% (3%) at an energy of 30 (100) GeV, and an angular^ 
resolution of about 0.01 radians. The energy response 141 
of the calorimeter to photons is calibrated using elec- 142 
trons from Z boson decays. S ince electrons and photons^ 
interact differently in the detector material before the 144 
calorimeter, additional energy corrections as a function 145 
of pj, are derived using a detailed GEANT-based simu- 146 
lation of the DO detector response. These corrections are i47 
largest, sa 2%, at photon energies of about 30 GeV. 148 

The data used in this analysis are collected using ai 49 
combination of triggers requiring a cluster of energy in™ 
the EM calorimeter with loose shower shape require-151 
mcnts. The trigger efficiency is «96% for photon candi-™ 
dates with p T «30 GeV and w 100% for p T > 40 GeV. ™ 

Offline event selection requires a reconstructed pp in-™ 
teraction vertex Q within 60 cm of the center of the™ 
detector along the beam axis. The missing transverse™ 



momentum in the event is required to be less than 0.7p T 
to suppress the background contribution from W —t ev 
decays. These requirements are highly efficient (> 98%) 
for signal events. 

The photon selection criteria in the current measure- 
ment are identical to those used in Ref. The pho- 
ton selection efficiency and acceptance are calculated 
using samples of 7 + c-jet events, generated using the 
SHERPA [lj| and PYTHIA [ll[ event generators. The sam- 
ples are processed through a GEANT-based [|| simulation 
of the DO detector response, followed by reconstruction 
usingthe same algorithms as applied to data. As in 
Ref. [5| , in the efficiency and acceptance calculations the 
photon is required to be isolated at the particle level by 
E™° = £5?* (0.4) -E T < 2.5 GeV, where £^ ot (0.4) is the 
total transverse energy of particles within a cone of ra- 
dius 1Z = \J (Ai]) 2 + (A0) 2 =0.4 centered on the photon 
and E T is the photon transverse energy. The particle 
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level includes all stable particles as defined in Ref. 
The photon acceptance varies within (82 — 90)% with a 
relative systematic uncertainty of (2 — 5)%, while the effi- 
ciency to pass photon identification criteria is (68 — 85)% 
with 3% systematic uncertainty. 

At least one jet with p^ > 15 GeV and |y jot | < 1.5 
must be reconstructed in each event. Jets are recon- 
structed using the DO Run II algorithm [l3| with a cone 
radius of TZ = 0.5. The jet acceptance with respect to 
the p£ and |y jGt | varies between 91% and 100% in differ- 
ent pip bins. Uncertainties on the acceptance due to the 
jet energy scale, jet energy resolution, and the difference 
between results obtained with SHERPA and pythia are 
in the range of (1 — 4)%. A set of criteria is imposed to 
have sufficient information to classify the jet as a heavy- 
flavor candidate: the jet is required to have at least two 
associated tracks with pr > 0.5 GeV with at least one hit 
in the SMT, and at least one of these tracks must have 
Pt > 1.0 GeV. These criteria have an efficiency of about 
90%. 

To enrich the sample with heavy-flavor jets, a neural 
net based 6-tagging algorithm (6-NN) is applied. It ex- 
ploits the longer lifetimes of 6-flavored hadrons in com- 
parison to their lighter counterparts, after the rejection 
of long-lived K° and A decays 



14|. The inputs to the 



6-NN combine information from the impact parameter 
of displaced tracks and the topological properties of sec- 
ondary vertices reconstructed in the jet to provide a con- 
tinuous output value that tends towards one for b jets 
and zero for light-quark jets. Events are required to con- 
tain at least one jet satisfying &-NN output > 0.7. This 
6-tagging selection suppresses light jets to less than 5% 
of the heavy-flavor enhanced sample. The efficiency for 
b and c jets to satisfy the ^-tagging requirements in the 
simulation is scaled by the data-to-Monte Carlo (MC) 
correction factors parametrized as a function of jet pt 
and r\ 14|. Depending on p T , the selection efficiency 
for this requirement is (8 — 10)% for c-jets with relative 
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systematic uncertainties of (6 — 23)%, caused by uncer- 
tainty on the data-to-MC correction factors. The maxi- 
mum difference between the efficiencies for c-jets arising 
from the Compton-like and annihilation subprocesses is 
about 10%. 

The relative rate of remaining light jets ("light/all") 
in the sample after the final selection is estimated us- 
ing SHERPA and pythia 7+jet events, taking into ac- 
count the data-to-MC correction factors as described in 
Ref. The light jet rates predicted by pythia and 

SHERPA agree within 5%. The central predictions are 
taken from SHERPA, which agrees with measured 7+jct 
15, Iff and 7 + fc-jet Q cross sections within (10 — 25)%. 
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After application of all selection requirements, 130,875 
events remain. We estimate the photon purity using an 
artificial neural network discriminant [H. The distribu- 
tion of the output of this discriminant (Onn) is fitted to a 
linear combination of templates for photons and jets ob- 
tained from simulated 7+jet and dijet samples, respec- 
tively. An independent fit is performed in each pip bin. 
It yields photon purities between 62% and 99%, which 
are close to those obtained in Ref. Their systematic 
uncertainties are of a comparable magnitude, (5-9)%. 

The invariant mass of all charged particles associated 
with a displaced secondary vertex in a jet, Msv, is a pow- 
erful variable to discriminate c from b jets. Since the Msv 
templates for light and c-jets after application of tight b- 213 
tagging requirements are quite close to each other, we 214 
first subtract the remaining small fraction (1 - 5%) of 215 
light jets from the data. Then the c-jet fraction is deter- 216 
mined by fitting Msv templates for c and b jets to the 217 
(7+heavy flavor jet) data. Jets from b quarks contain see- 218 
ondary vertices that have in general larger values of Msv 219 
as compared to c jets and the region beyond Msv > 2.0 220 
GeV is strongly dominated by b jets. The templates for 6 221 
and c jets are obtained from pythia samples of 7 + 6-jet 222 
and 7 + c-jet events, respectively, and are consistent with 223 
the templates generated using SHERPA. The templates 224 
for jets arising from the Compton-like and annihilation 22 ^ 
subprocesses are also similar to each other. 226 

The result of a maximum likelihood fit to the Msv 227 
templates, normalized to the number of events in data, 228 
is shown in Fig. Q] for the 50 < pip < 60 GeV bin as an 22 ' 
example. Fits in the other pip bins are of similar quality. 230 
As shown in Fig. ^ the estimated c-jet fraction obtained 231 
from the fits in the final selected heavy-flavor sample af- 232 
ter subtraction of the light-jet component drops with in- 233 
creasing p T , on average, from about 52% to about 40%. 234 
The corresponding fit uncertainties range between (4- 23 5 
32)%, increasing towards higher p T , and are dominated 2 ^ 
by the limited data statistics. Since the fits are per- 23 ? 
formed independently in each pip bin, these uncertainties 233 
are uncorrelated from bin to bin. Additional systematica 
uncertainties are estimated by varying the relative rate ofiio 
light jets by ±50% and by considering the differences in24i 
the light jet predictions from SHERPA and PYTHIA event242 



> 

(D 
10 

2-4 

w 

c 

0) 

m 3 





* 

□ 
1 




n(7J 1 - ft 7 fh" 1 


- 




I 

T 


• Data 






b-jet 






c-jet 


- 




-a- b+c jets 








50 < p y < 60 GeV 






in 


1 IlS'lkm.i.a.tiUa.j,, 



"0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
Secondary Vertex Mass (GeV) 

FIG. 1: (Color online) Distribution of secondary vertex mass 
after all selection criteria for a representative bin of 50 < < 
60 GeV. The expected contribution from the light-jet compo- 
nent has been subtracted from the data. The distributions 
for the fo-jet and c-jet templates (with statistical uncertain- 
ties) are shown normalized to their respective fitted fractions. 



generators. These two sources lead to uncertainties on 
the c-jet fraction of about (5-9)% and 6%, respectively. 

Systematic uncertainty on the measured cross sections 
due to the 6-NN selection is estimated by performing the 
measurement with looser 6-NN selections: requiring 6-NN 
output > 0.3 or > 0.5 instead of 0.7. In both cases, this 
significantly increases the light-jet rate and also changes 
the c- and 6-jet fractions, resulting in a variation of the 
7 + c-jet cross section of < 7%. This variation is taken 
as a systematic uncertainty on the cross section. 

The data, corrected for photon and jet acceptance, 
reconstruction efficiencies and the admixture of back- 



ground events, are presented at the particle level [12j for 
comparison with predictions by unfolding the data for 
effects of detector resolution. 

The differential cross sections of 7 + c-jet production 
are extracted in nine bins of p T . They are listed in Ta- 
ble U and are shown in Fig. [3] The data points are plotted 
at the values of pip for which the value of a smooth func- 
tion describing the dependence of the cross section on pip 
equals the averaged cross section in the bin [l?} . 

The statistical uncertainty of the results ranges from 
2% in the first p T bin to 11% in the last pip bin. The 
total systematic uncertainty varies between 14% and 42% 
across these bins. The main sources of uncertainty at low 



pip are due to the photon purity (up to 



the c-jet 



fraction (10 - 33%), and the luminosity (6%) 0. The 
total systematic uncertainties (5 S yst) and the bin-to-bin 
uncorrelated components (5"^) are shown in Table |U 
Next-to-leading order (NLO) perturbative QCD pre- 
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TABLE I: The 7 + c-jet production cross sections da/dp T in bins of pip for \y 1 \ < 1.0, > 15 GeV and \rf et \ < 1.5 together 
with statistical uncertainties (<5 s tat), total systematic uncertainties (<5 sy st), and the uncorrelated component of S syst (S^yst)- The 
column Stot shows total experimental uncertainty obtained by adding 5 sta t and <5 sy st in quadrature. The last four columns show 
theoretical predictions obtained within NLO QCD, /cr-factorization, and by the pythia and SHERPA event generators. 
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FIG. 2: The c-jet fraction in data after subtraction of light- 
jet background as a function of plj, derived from the template 
fit to the heavy quark jet data sample after applying all se- 
lections. The error bars include statistical and systematical 
uncertainties. Binning is the same as given in Table [U 



dictions of order O (a^ ) [H, EH , with the renormalization 
scale fin, factorization scale jip, and fragmentation scale 
Hf all set to pj,, are given in Table HI The uncertainty 
from the scale choice is estimated through a simultane- 
ous variation of all three scales by a factor of two, i.e., 
for fJ-R.F.f = 0.5pJ and 2p7p, and is found to be similar 
to those for 7 + 6-jet predictions (5 — 30)%, being larger 
at higher p\ The NLO predictions utilize CTEQ6.6M258 
PDFs [l9[ and are corrected for non-perturbative effectS259 
of parton-to-hadron fragmentation and multiple parton26o 
interactions. The latter are evaluated using SHERPA261 
and pythia MC samples generated using their default262 
settings [HI EH. The overall corrections vary withiii263 
0.90 — 0.95 with an uncertainty of < 2% assigned to ac-264 
count for the difference between the two MC generators. 265 
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FIG. 3: (Color online) The 7 + c-jet differential production 
cross sections as a function of p T . The uncertainties on the 
data points include statistical and systematic contributions 
added in quadrature. The horizontal error bars show the p T 
bins. The measurements are compared to the NLO QCD cal- 
culations [l|,|lil using CTEQ6.6M PDFs [3 (solid line). The 
predictions from SHERPA [l^, pythia [ill] and kr factoriza- 
tion approach [2^, [2l[ are shown by the dash-dotted, dotted 
and dashed lines, respectively. 



The predictions based on the /cT-factorization ap- 
proach 2^, 21 1 and unintegrated parton distributions [22[ 
are also given in Table HI The resummation of gluon di- 
agrams with gluon transverse momentum (&t) above a 
scale \x of order 1 GeV, leads to a broadening of the 
photon transverse momentum distribution in this ap- 
proach . The scale uncertainties on these predictions 
vary from about -28%/ + 31% at 30 < < 40 GeV to 
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FIG. 4: (Color online) The ratio of 7 + c-jet production cross3i3 
sections to NLO predictions for data and theoretical predic-314 
tions. The uncertainties on the data include both statistical 315 
(inner error bar) and total uncertainties (full error bar). Also 3l 
shown are the uncertainties on the theoretical QCD scales 
and the CTEQ6.6M PDFs. The ratio for intrinsic charm mod- 
els [2^ ] are presented, as well as the prediction s gi ven by 318 
fcx-factorization [HH]]], SHERPA [l(J and pythia [lj|. 319 



about +14%/ + 5% in the last p T bin. 322 
Table U also contains predictions from the pythia fill] 323 
event generator with the CTEQ6.1L PDF set. It includes324 
only 2 — > 2 matrix elements (ME) with gc — > 7c and325 
qq — > 7.g scatterings (defined at LO) followed by g — > cc™ 
splitting in the parton shower (PS). We also provide pre-327 
dictions by the SHERPA MC event generator [n| witli328 
the CTEQ6.6M PDF set [12]. Matching between the ME329 
partons and the PS jets follows the prescription given in 3 3o 
Ref. [l5|, with the matching scale taken to be 15 GcV.331 
Systematic uncertainties are estimated by varying thc332 
ME-PS matc hing scale by ±5 GcV around the chosen333 
central value 
ation. 



231 ] . resulting in a ±7% cross section vari- 



All theoretical predictions are obtained using the pho-336 
ton isolation requirement of Ep° < 2.5 GcV. The pre- 337 
dictions are compared to data in Fig. [3] as a function of338 
pTp. The ratios of data over the NLO QCD calculations339 
and of the various theoretical predictions to the NLO340 
QCD calculations are presented in Fig. The NLO341 
predictions with CTEQ6.6M agree with MSTW2008 0342 
and ABKm09nlo [25} within 10%. Parameterizations343 
for models containing intrinsic charm (IC) have been in-344 
eluded in CTEQ6.6c j26(. Here we consider the BHPS IC345 
model [13, 28 1, based on the Fock space picture of the346 
nucleon structure (29|, in which intrinsic charm appears34? 



mainly at large momentum fractions x, and the sca-likc 
model in which the charm PDF is sea-like, similar to 
that of the light-flavor sea quarks. The NLO QCD pre- 
dictions based on these intrinsic charm models are nor- 
malized to the standard CTEQ predictions and are also 
shown in Fig. 0] Both non-perturbative intrinsic charm 
models predict a higher 7 + c-jet cross section. In the 
case of the BHPS model, the ratio grows with p T , while 
an opposite trend is exhibited by the sea-like model. 

The measured cross sections are in agreement with the 
NLO QCD predictions within theoretical and experimen- 
tal uncertainties in the region of 30 < p T < 70 GeV, 
but show systematic disagreement for larger p T . The 
cross section slope in data differs significantly from the 
NLO QCD prediction. The results suggest a need for 
higher-order perturbative QCD corrections in the large 
Pp region, which is dominated by the annihilation process 
qq — > 7(7 (with g — > cc), and resummation of diagrams 
with additional gluon radiation. In addition, the underes- 
timation of the rates for diagrams with g — > cc splittings 
may result in lower theoretical predictions of cross sec- 
tions as suggested by LEP 0, LHCb [HI and ATLAS 
[32^ results. The prediction from the Aix-factorization ap- 
proach is in better agreement with data at p T > 120 GeV. 
However, it underestimates the cross section in the low 
and intermediate p T region. The 7 + c-jet cross section 
as predicted by SHERPA becomes higher than the NLO 
QCD prediction at large p^,, but is still lower than the 
measured values. It has been suggested that combining 
NLO parton-level calculations for the ME with PS pre- 
dictions [33| will improve the description of the data [34| . 

In addition to measuring the 7 + c-jet cross-section, 
we also obtain results for the 7 + 6-jet cross section us- 
ing the new tight 6-NN selection. The values of the 
obtained 7 + fe-jct cross section agree within 10% (i.e. 
within uncertainties) with the published results [5] ob- 
tained with a looser 6-NN selection. We use them to 
calculate the ratio a(j + c)/a{^ + b) in bins of p T . In 
this ratio, many experimental systematic uncertainties 
cancel. Also, theory predictions of the ratio are less sen- 
sitive to the scale uncertainties, and effects from missing 
higher-order terms that impact the normalizations of the 
cross sections. The remaining uncertainties are caused by 
largely (65 — 67%) correlated uncertainties coming from 
the fitting of c-jet and 6-jet Msv templates to data, and 
by other uncertainties on the c-jet fractions discussed 
above. The systematic uncertainties on the ratio vary 
within (6 — 26)%, being largest at high p T . Theoreti- 
cal scale uncertainties, estimated by varying scales by a 
factor of two (to Hr,fj = 0-5p T and 2p T ) in the same 
way for (7(7 + c) and (7(7 + b) predictions, are also signif- 
icantly reduced. Specifically, residual scale uncertainties 
are typically < 10% for the fex-factorization approach 
and < 4% for NLO QCD, which indicates a much smaller 
dependence of the ratio on the higher-order corrections. 
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TABLE II: The a(-y + c)/a(-y + b) cross section ratio in bins of pj for |y 7 | < 1.0, p£* > 15 GeV and |r; jot | < 1.5 together 
with statistical uncertainties (<5 s tat), total systematic uncertainties (<5 sy st), and the uncorrelated component of S syst (^syst)- The 
column 5 to t shows total experimental uncertainty obtained by adding 5 sta t and <5 sy st in quadrature. The last four columns show 
theoretical predictions obtained using NLO QCD, /ct -factorization, pythia and SHERPA event generators. 
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FIG. 5: (Color online) The ratio of 7 + c-jet and 7 + 6-jet 
production cross sections for data together with theoretical 378 
predictions as a function of p^,. The uncertainties on the data 379 
include both statistical (inner error bar) and total uncertain- 380 
ties (full error bar). Predictions given by fcx-factorization 
[U EI, SHERPA [!(} and pythia [L^ are also shown. The 381 
pythia predictions with a contribution from the annihilation 382 
process increased by a factor of 1.7 are shown as well. The 383 
predictions for intrinsic charm models [2(| are also presented. 384 



Experimental results as well as theoretical predictions for 387 
the ratios are presented in Table [XXJ 3 as 
Figure [5] shows the measured ratio <7 (7 + ^/(7(7 + 6) as389 
a function of pj, and a comparison with various predic-390 
tions. There is good agreement with NLO QCD, SHERPA391 
and pythia predictions in the region 30 < pip < 7O392 
GeV, while /ci-factorization predicts smaller ratios than 3 9 3 
observed in data. At higher p T , data show systemati-394 



cally higher ratios than NLO QCD, SHERPA and pythia 
predictions, while fex-factorization starts agreeing with 
data within uncertainties. We also show NLO predic- 
tions with the BHPS 



27 



28j and sea-like IC models |2fJ 
used to predict 7 + c-jet cross section, while standard 
CTEQ6.6M is used to predict the 7 + &-jet cross section. 
The BHPS model agrees with data at p T > 80 GeV, 
while the sea-like model is significantly beyond the range 
of data points. BHPS model would better describe the 
ratio to data with a small shift in normalization. As 
with the 7 + c-jet measurement, the 17(7 + c)/er(7 + b) 
ratio can also be better described by larger g — > cc rates 
than those used in the current NLO QCD, SHERPA and 
pythia predictions. To test this, we have increased the 
rate of the annihilation process (where c jet is always pro- 
duced due to g — > cc splitting) in the PYTHIA predictions. 
The best description of data is achieved by increasing the 
rates by a factor of 1.7 with x 2 /ndf ~ 0.7 (compared to 
X 2 /ndf = 4.1 if such a factor is unity). However, accord- 
ing to our estimates using the signal events simulated 
with SHERPA, there are also about (10-35)% (higher for 
larger p T ) events with two c-jets. Assuming that one jet 
is coming from gluon initial state radiation followed by 
g — > cc splitting, the required overall correction factor 
would be smaller by about (8-24)%. 

In conclusion, we have measured the differential cross 
section of 7 + c-jet production as a function of p^ at 
the Fermilab Tevatron pp collider. Our results cover the 
kinematic range 30 < pi < 300 GeV, p£ > 15 GeV, 
I y 1 1 < 1.0, and \y iCt \ < 1.5. In the same kinematic re- 
gion, and in the same p^ T bins, we have measured the 
(7(7 + c)/a(-f + b) cross section ratio. None of the theo- 
retical predictions considered give good description of the 
data in all p^, bins. Such a description might be achieved 
by including higher-order corrections into the QCD pre- 
dictions, while at > 80 GeV the observed difference 
from data may also be caused by an underestimated con- 
tribution from gluon splitting g — > cc 33"32| in the anni- 
hilation process or by contribution from intrinsic charm. 
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The presented results can be used for further develop-436 
ment of theoretical models to understand production of 437 
high energy photons in association with heavy flavor jets. 438 
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